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Abstract 

Here, creation of the universe is obtained only from gravity sector. 
The dynamical universe begins with two basic ingredients (i)vacuum 
energy, also called dark energy (as vacuum energy is not observed) and 
(ii) background radiation. These two are obtained through one-loop 
renormalization of riccion. Solutions of renormalization group equa- 
tions yield initial value of vacuum energy with density p\ ew = 10 6 GeV 4 
as well as show a phase transition at the electroweak scale M ew . As a re- 
sult of phase transition, energy is released (in the form of background 
radiation) heating the universe upto temperature T ew = 78.5GeV = 
9.1 x 10 14 i^ initially. In the proposed cosmology, it is found that not 
only current universe expands with acceleration, but it undergoes ac- 
celerated expansion from the beginning itself. It is demonstrated that 
dark energy decays to dark matter as well as ratio of dark matter 

density and dark energy density remains less than unity upto a long 

l 
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time in future universe also, providing a solution to cosmic coincidence 
problem. Future course of the universe is also discussed here. It is 
shown how entropy of the universe grows upto 10 87 in the present 
universe. Moreover, particle creation, primordial nucleosynthesis and 
structure formation in the late universe is discussed for the proposed 
model. Thus, investigations, here, present a fresh look to cosmology 
consistent with current observational evidences as well as provide so- 
lution to some important problems. 
PACS nos.98.80.Cq, 04.62. +v, 04.50. +h, 95.35.+d. 
Key Words: Origin of the universe, Dark Energy and Dark Matter, 
Higher-dimensional higher-derivative gravity, Quantum Field Theory 
in curved spaces, One-loop renormalization. 



1. Introduction 

Experimental probes, like luminosity measurements of high redshift 
supernova [1], anisotropies in cosmic microwave background [2] and ob- 
servational gravitational clustering [3], strongly indicate late time ac- 
celerated expansion of the universe. Theoretically, accelerated expan- 
sion of the universe can be obtained either by modifying left hand side 
of Einstein's equations, like introduction of C-field in steady— state 
theory adhered to the Perfect Cosmological Principle [4] or by using 
energy momentum tensor , having dominance of exotic matter, with 
negative pressure violating the strong energy condition. This kind of 
matter is known as dark energy (vacuum energy), which has drawn 
much attention of cosmologists today. Past few years have witnessed 
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concerted efforts to propose different dark energy models. Important 
efforts, in this direction, are scalar field models like (i) quintessence 
[5], (ii) k-essence [6], (iii) tachyon scalar fields [7, 8, 9, 10, 11, 12] and 
models based on quantum particle production, Chaplygin gas [13, 14]. 
In these models, a scalar field acts as a source of are dark energy and 
plays crucial role in the dynamical universe. But these models have 
no answer to the question "Where are these scalars coming from?" It 
is like Higg's fields in GUTs as well as inflaton in inflationary models 
of the early universe. In the cosmology probed here, no scalar field is 
required to be incorporated, from outside the theory, to discuss cosmic 
dark energy. Following arguments in [14, 15], here also, vacuum energy 
is recognized as DE. 

All these dark energy models (mentioned above) try to explain rolling 
down of DE density from a very high value in the early universe to 
an extremely small value in the current universe as it is suggested 
by astronomical observations. But its initial value, in the early uni- 
verse, is different in different contexts. For example, it is ~ 10 76 GeV 4 
at Planck scale, ~ 10 60 GeV 4 at GUT phase transition and quantum 
chromodynamics yields its value ~ lCT 3 GeV 4 . In what follows, one- 
loop renormalization of riccion ( a particle representing physical role 
of the Ricci scalar explained below) contributes initial value equal to 
10 6 GeV 4 for slowly varying time-dependent DE density to the dynami- 
cal universe beginning with a phase transition at the electro-weak scale 
M ew = lOOGeV. 

Work starts effectively from a general case of (4 + D) -dimensional 
space-time and, later on, it is shown that D — > 6 making effective 
dimension of the space-time equal to 10. The observable universe is 
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a 4-dimensional hypersurface of the higher-dimensional world. So far 
sharpest experiments could probe gravity upto 0.1 mm. and, in the 
matter sector, probe could be possible upto electroweak scale M ew = 
lOOGeV, having wavelength as short as ~ 1.97 x lCT 16 cm. It means 
that it is difficult to realize 4-dimensional gravity for length scales less 
than 0.1 mm. So, it is reasonable to think higher-dimensional gravity 
for scales smaller than this scale. This idea is parallel to the brane- 
world gravity, where it is assumed that gravity is stronger in higher- 
dimensional space-time, called bulk and only a small part of it is relized 
in the observable universe [17]. 

It has been noted by many physicists that the Ricci scalar R behaves 
like a physical field also, in addition to its geometrical nature, if gravi- 
tational action contains higher-derivative terms [18, 19, 20, 21, 22, 23, 
24, 25]. In 1980, Starobinsky suggested that if sign of the R 2 term in the 
higher-derivative gravitational action is chosen properly, one could ob- 
tain only one scalar particle with positive energy and positive squared 
mass. He called it as "scalaron" [18]. In papers [18, 19], gravitational 
constant G is either taken equal to unity or lagrangian density is taken 
as j^q(R+ higher - derivative terms). As a result, (mass) 2 of the Ricci 
scalar, does not depend on the gravitational constant G. Realizing the 
important role of G in gravity, in papers [20, 21, 22, 23, 24, 25] as 
well as here, lagrangian density is taken as ( w 1 7tG R+ higher - deriva- 
tive terms) leading to a drastic change where (mass) 2 depends on G 
also. In the following section as well as earlier works , it is shown that 
physical aspect of the Ricci scalar is given by a scalar field R = r/R {rj 
is a parameter having length dimension), called riccion with (mass) 2 
depending on gravitational constant G and other coupling constants in 
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the action [20, 21, 22, 23, 24, 25]. It is different from scalar on in two 
ways (i)mass dimension of riccion is one like other scalar fields ,such 
as quintesence, inflaton and Higg's scalar, whereas mass dimension of 
scalaron is two and (ii) (mass) 2 for scalaron does not depend on G, 
whereas ,for a riccion, it depends on G. 

The present paper begins with higher-derivative gravitational action 
in higher-dimensional space-time with topology M 4 ® S D , where S D is 
a D-dimensional sphere being the hidden extra-dimensional compact 
space. Higher-derivative gravity faces ghost problem, which can be 
avoided by taking coupling constants in the action properly. 

The distance function for (4 + D) -dimensional space-time is defined 

as 

dS 2 = g^dx»dx u - l 2 dtt 2 (1.1a) 

with 

dtt 2 = d9l + s%n 2 e x del + ■■■ + sin 2 6 1 ■ ■ ■ sin 2 e {D - X) de 2 D . (1.16) 

Here ^(/i, v — 0, 1, 2, 3) are components of the metric tensor in M 4 , 
/ is radius of the sphere which is independent of coordinates x^ and 
< 0i, 2 , • • • , (D -i) < 7T and < D < 2vr. 

It is important to mention here that riccion is different from the 
scalar mode of graviton. which is highlighted in Appendix A. 

Thus, in the proposed cosmological scenario, our dynamical universe 
begins with the initial dark energy density (vacuum energy density)pA ow = 
1.16 x 10 7 GeV 4 and background radiation with temperature T ew = 
33.5GeV = 3.89 x lO 14 ^, contributed by riccion at M cw . The back- 
ground radiation is caused by a phase transition at M cw . This event 
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is recognized as big — bang. As usual, dynamical universe grows in the 
post big — bang era. In the present theory, riccion has effective role in 
the pre big — bang era, but it remains passive in the post big — bang 
times. 

The paper is organized as follows. Section 2 demonstrates derivation 
of riccion equation and its action. One-loop renormalization of riccion 
is done in sestion 3. Renormalization group equations are solved in 
section 4. This is an important section as it provides initial values of 
dark energy density pa and temperature T of the observable universe. 
Moreover, in this section, dimension of hidden space is derived to be 6. 
In section 5, derivation of equation of state for dark energy (vacuum 
energy) is derived and phase transition at M ew is discussed. It is found 
that dark energy decreases with expansion of the universe. So, it is 
natural to think for decay of dark energy to dark matter, which is dis- 
cussed in section 6. This section is very important from cosmological 
point of view, as many important results are derived here. It is demon- 
strated that dark energy decays to hot dark matter (HDM) till matter 
remains in thermal equilibrium with radiation, but when temperature 
falls down the decoupling temperature, production of HDM decreases 
and creation of cold dark matter (CDM)increases. Dark matter den- 
sity is found less than dark energy density from the epoch of big — bang 
upto the time 10.6to ( is the present age of the universe ). This re- 
sult provides a solution to cosmic coincidence problem, which raises 
the question " Why does dark energy dominates over dark matter 
only recently?." First time, this question was posed by P.J.Steihardt 
[26]. In some earlier works also, this problem was adrressed and so- 
lutions were suggested, taking coupled system of quintessence scalar 
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fields and matter [27, 29, 30, 31]. As mentioned above, contrary to 
earlier attempts, such scalar fields are not required here. Moreover, it 
is found that the universe undergoes an accelerated expansion from the 
beginning itself upto late future universe. Future course of dynamics 
is also discussed here. In this section, temperature and entropy of the 
universe are discussed and it is demonstrated how entropy of the uni- 
verse grows upto 10 87 upto the current epoch. In section 7, creation of 
spinless and spin- 1/2 particles , primordial nucleosynthesis as well as 
structure formation in the late universe are discussed. The last section 
summarizes results. 

Natural units, defined as % = /i = c = 1 ( where Kb is Boltzman's 
constant, H is Planck's constant divided by 2n and c is the speed of 
light) are used here with GeV is used as a fundamental unit such that 
lGeV = 1.16 x 10 13 K = 1.78 x l(T 24 #m , lGeV" 1 = 1.97 x l(T 14 cm = 
6.58 x 1(T 25 sec. 

2.Riccions from (4+D)-dimensional geometry of the space-time 

Theory begins with the gravitational action 

d 4 xd D yy/-g^ +D ) — ^(4+d)#(4+d)+7(4+d)(#(4+d) 



6P + 3) \ 

" (£) _ 2) ( D + 4 ) K (D+4); 



(2.1a) 



where G(4+d) = M _ ( 2+D ) (M being the mass scale ),oc(a+d) = otVp 1 , 7(4+d) = 
3! (d- 2) ^0 1 ' an d Rd = D ^i~^ ■ Here Vd, being the volume of S D ,is 
given as 

27r (D+l)/2 

V D = — —l D . (2.1b) 

r(D + i)/2 v ; 
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gu+D) is determinant of the metric tensor qmn{M, N = 0, 1, 2, • • • , (3 + 
D) and R{a+d) — R + Rd- a is a dimensionless coupling constant, R is 
the Ricci scalar in M 4 and Gn = G^ + d)/Vd- 

Invariance of S under transformations ()mn — > 9mn + &9MN yields 
[24, 25] 



M( 2+D ) 1 (1) (2) 

— ^ (Rmn — -gMNR(4+D)) + ot(4+D)H MN +'y(4+D)H MN = 0, (2.2a) 
where 



HltfN ~ 2-R;MAf ~ < 2gMN a (4+D)R(4+D) ~ ^fi'MAf-R^+D) + ^R(4+D)RmN, 



(2.26) 



and 
H 



(2) oo2 o. rn n2 6(D + 3) ( 1 



MiV - 3R. MN -3gMN D (4+D)R(4 + D)- (_o _ 2) { _ 2 fl,MArD(4+jD) " R(4+D) 

+ 2n(4 +D )_R(£) + 4)_R MA r + R 2 MN } — ydMNRfl+D) + ^R 2 (A+D)^MN (2.2c) 

with semi-colon (;) denoting curved space covariant derivative and 
Trace of these field equations is obtained as 



R(4+D) + «(4+D) [2{U + 3J U(4+D)R(4+D) + ^R{4+D)\ 



32tt 

+i 7(4+D) (D-2)i?f 4+D) = 0. (2.3) 
In the space-time described by the distance function defined in eq.(l.l), 



□ (1+D) fi (1+D) = M = 7 LA( V ^ 9 ,^) fl , ( 2 . 4) 
using the definition of R(4+d) given in eq.(2.1). 
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Connecting eqs.(2.3)- (2.4) as well as using R^+d), &(a+d) an d 7(4+d) 
from eq.(2.1), one obtains in M 4 

_ [ (J + 2)^»V D ](fi + flo)+a|2(fl + 3)Dfl 

+\d{R + Rd) 2 ] + ^( R + r d? = 0, (2.5) 
which is re-written as 

[□ + -iR + m 2 + -rfR^R + r]-^ = 0, (2.6) 
2 3! 

where 



£ = 



m 



A 



2(D + 3) + ^ 
(D + 2)\M^V D 
16tt 



2(D + 3) 2 ' D 



4( J D + 3)a' 

(.D + 2)AM( 2+D )\/ D 



16tt 



+ 



DRl 



4(L> + 3) 6' D J 



(2.7a, 6, c, d) 

where a > to avoid the ghost problem. 

A scalar field, representing a spinless particle, has unit mass di- 
mension in existing theories. R, being combination of second order 
derivative as well as squares of first order derivative of metric tensor 
components with respect to space-time coordinates, has mass dimen- 
sion 2. So, to have mass dimension like other scalar fields, eq.(2.6) is 
multiplied by rj and r\R is recognized as R. As a result, this equation 
looks like 
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[n + ^R + m 2 + ^-R 2 ]R + ^ = 0. (2.8) 

For scalar on, n is dimensionless [18]. 

The above analyses show that, on taking trace of field equations (2.2) 
and compactifying the space-time M 4 ® S D to M 4 , only one degree 
of freedom is obtained spontaneously, which is the scalar mode R [16, 
18-23] . It is unlike gravitons having 5 degrees of freedom including one 
scalar. In Appendix A, it is explained that scalar mode of graviton is 
different from riccion. 

If R is a basic physical field, there should be an action yielding 
eq.(2.8) for invariance of S^, under transformations R — > R + SR. 

In what follows, is obtained. If such an action exists, one can 
write 

5S k = - J d 4 x^g5R[(a + ^R + m 2 + ^R 2 )R + &] (2.9a) 

which yields eq.(2.8) if 5Sp, = under transformations R — > R + 5R. 
Eq.(2.9a) is re- written as 

SS A = J d 4 x^^Rd^(5R) - [^RR 2 + m 2 R + ^R 3 + &)6R 

= J d 4 x5\y=g~\^d> x Rd lx R - (^RR 2 + \m 2 R 2 + ^R 4 + &r) 

(2.96) 

R, R and d 4 x^/—g are invariant under co-ordinate transformations. 
So, R(x) = R(X),R(x) = R(X) and d 4 x^g = d 4 X, where X\i = 
0, 1, 2, 3) are local and x l (i = 0, 1, 2, 3) are global coordinates. More- 
over, 



U - g ii 92 | 1 g mnd9mn gij d | dg* d d 2 



dx l dxi 2 dx l dxi dx l dxi dX l dXi 
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in a locally inertial co-ordinate system, where gij = rjij (components 
of Minkowskian metric) and g l \ = (comma (,) stands for partial 
derivative). Thus, in a locally inertial co-ordinate system, 

6S A = J d 4 x8{V=g-[\d»Rd,R - (^RR 2 + \rn 2 R 2 + ±fr + tfi?)] } 
= J d'X 6{y/=j[±9'Rd ll R - (^RR 2 + \m 2 R 2 + ^R 4 + 
= 5 J d A X {y/^jj^&RduR - (^RR 2 + \m 2 R 2 + ^R 4 + &R) 

Employing principles of covariance and equivalence as well as eq.(2.9b) 
which implies that 

S k = j d 4 x{^9~[^Rd,R - QjZRR 2 + \m 2 R 2 + ^R 4 + w)] }. 

(2.10) 

It is important to mention here that R is different from other scalar 
fields due to dependence of (mass) 2 on the gravitational constant, di- 
mensionality of the space-time and the coupling constant a, given by 
the eq.(2.7b). Moreover, it emerges from geometry of the space-time. 



3. One-loop quantum correction and renormalization of riccion 

The S k with the lagrangian density, given by eq.(2.10), can be ex- 
panded around the classical minimum R in powers of quantum fluc- 
tuation R g = R — R as 



'R 



Sf + + sf + 



R 



R 



R 
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where 
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/ 



sf 



and 



1 

A 



d A x^-g [-g^dpRodvRo - (3^0 + 2 m ^o + + ^o) 



d 4 Xy/^R q [D + -£R + m 2 + -R 2 }R 



*«=0 

as usual, because this term contains the classical equation. 

The effective action is expanded in powers of h (with h — 1) as 

r(R) = s R + r {1) + r' 

with one-loop correction given as [32, 33] 



= -lnDet{D/n 2 ), 



(3.1a) 



where 



D 



6 2 S R 



5R 2 



R=Ro 



(3.16) 



and T' is a term for higher-loop quantum corrections. In eq.(3.1), ji is 
a mass parameter to keep dimensionless. 

To evaluate the operator regularization method [34] is used upto 
adiabatic order 4. As potentially divergent terms are expected upto this 
order only. In a 4-dim. theory, one-loop correction is obtained as 



r« = 



_2\-l ^ 



167r T" 
as 



d 4 x v / -g(x)(^ T 



M 4 



M 2 



1 1 



£ )R + 



-(- 

L6 V5 



1 



1) V6 



(a-2)(a-l) 
1 



s=0 
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(3.2a) 

where 

M 2 = m 2 + (X/2)R 2 . (3.26) 

Here it is important to note that matter as well as geometrical both 
aspects of the Ricci scalar are used in eq.(3.2). The matter aspect is 
manifested by R and the geometrical aspect by R, Ricci tensor compo- 
nents R^ v and curvature tensor components R^uap as it is mentioned 
above also. 

After some manipulations, the lagrangian density in is obtained 

as 

Lrm = ( i«o-[(^ + (V2)ja'{|-^( ^±^g )} 

180 ' 2\6 2 1 II 



(3.3) 

Now the renormalized form of lagrangian density can be written as 



Lren = l -g» v d^R^R* - ^R 3 - \m 2 R\ - ^R* - m + A 
+e 4 DR + L r( i) + L ct 



(3.4a) 
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with bare coupling constants A« = (m 2 , A, A, £, eo, ei, e2, 63, £4), I^ 1 
given by eq.(3.3) and L ct given as 



-i^^o - 7^™ 2 ^o - ijj^o - 5 $Ro + SA + 5e R + \^R 2 



+5e 2 R^R^ + SesR^R^ + 5e A UR. 



(3.46) 



In eq.(3.4b), 5\ = (5m 2 , 5$, 5\, 5A, <5£, 5e , Sei, 5e 2 , Se 3 , o~e 4 ) are counter- 
terms, which are calculated using the following renormalization condi- 
tions [35, 36] 



A 
A = 

= 

m 2 = 



£4 



= L 



& 



<rcn\ Ru= R (0)0>R= 

=-R(0)l,fl=0 

rt() 

flo=fl (0) i,fl=0 



a" 

~dR 
d 2 



,4 -^ ren 



v rcn 



dRl 



^ren 



-77 



03 



_R =0,_R=0 

L. 



R =R (0)2 ,R=0 



d_ 
dR 

d 2 



L 



ren 



R =0,i?=0 



dR? R =0,R=R 5 
d(R^R^) rcn R =0,R=R 6 

9 L 

9 , 



=#7 



9(Di2)" 



y ren 



_R o =0,R=i?g 

(3.5a, 6, c,d,e, f,g,h,i,j) 
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As R = nR, so when R = 0, R(o)o — R(o)i — R(o)2 = and R$ = 
R G = R 7 = R 8 = when Rq = 0. 

4.Renormalization group equations and their solutions 

Mass scale dependence of coupling constants are obtained by solving 
renormaliztion group equations 



d\ 
dr 



(4-1) 



where r = \ln(M^ / fj, 2 ) and (5^ are one-loop f3 functions for different 
coupling constants. Here \i is the mass scale parameter and M ew is 
the cut-off scale such that /i > M cw . As experimental probes could be 
possible upto M ew , so this scale is used as a cut-off mass scale. 

Px i in eq.(4.1) are calculated using counter-terms yielded by eqs.(3.4) 
and (3.5) and putting /i^Aj = for bare coupling constants in equa- 
tions 



(4.2) 



Using /5-functions for different coupling constants, given by eqs.(4.2), 
solutions of differential equations (4.1) are derived, as 



A 
A 



m 



A + mew 

1 'cw \ 



= A 



" w ' 2A 
1 - 



3A ew r\i/3 



/ _ £AewT\ 

V 8tt 2 J 



ew 

3A nw ri - 1 



8tt 2 

$ ew (constant) 
3A Rw r 



1 - 



-1/3 



6 V2 W 6J I 8tt 2 J 



(4.3a, b, c, d, e) 
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for coupling constants of relevant terms to be used in further inves- 
tigations in the present paper. Here, Aj ow = Aj(r = 0) and r = 
at (j, — M ew according to defnition of r given above. 

These results show that as pL — > oo(r — > — oo), A — > and m 2 — > 
and |e -> I- 

Using these limits in eq.(2.7a), it is obtained that 



D = 6. 

Also,eqs.(2.1a) and (4.4) imply 



(4.4) 



Re, 



30 



So, from eq.(4.3a) 



A = Ar 



(4.5) 



(4.6) 



at /j, — M ew . The equation (4.3c) shows that d is independent of mass 
scale ji. So, rj^ lf d, being true for all /x, is obtained from eq.(2.7d) as 



30 



'/ » n 



m ew- 



5(1 + 60A C 



at /i — M ew . Here 



m 



A cw M c 8 w U 6ew , 10(l + 450A e 



2tt 



+ 



l ew 



with 



1 fi7T 3 / 6 
t 7 _ ±u/l l ew 



6cw 
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(4.6) 



(4.7) 



(4.8) 



which is derived connecting eqs.(2.7b,d), taking the arbitrary parame- 
ter 7] = l cw and putting D = 6. 

Thus rj~ lf d , given by eq.(4.6),is an imprint of extra six-dimensional 
compact manifold S 6 in the 4-dimensional universe. This term has the 
dimension of energy density. Moreover, it is not generated through 
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matter, but the geometry of extra-dimensional space. So, it is recog- 
nized as dark energy density pa 6W at fj, — M cw , given as 



30 



I 



2 

ew 



5(1 + 60A C 



PA CW = IT 1 * = ^ ml - ^tZ^HL , (4.9) 



2 
cw 



5. Equation of state for dark energy and Phase transition 
at the eletroweak scale 

Dark energy density Pa ow , obtained through renormalization of ric- 
cion at energy mass scale \x = M ew , is given by eq.(4.9). As riccion 
is emerging from geometry of the space-time, pa cw can be obtained as 
zero-point energy of riccion also, given as 

p A (t = 0) = p Aow = (2tt)- 3 / y/k^ + ml w A7Tk 2 dk, (5.1) 

Jo 

where m ew is the mass of riccion at the electroweak scale M ew . Eq.(5.1) 
shows that pa cw diverges as k — * oo. But pa cw == A cw (as given by 
eq.(4.7)), being finite, implies that the integral in eq.(5.1) should be 
regularized upto a certain cut-off mode k = k c . As so far experiments 
could be performed upto M ew only, so cut-off scale is taken as k c = M ew 
as above. Moreover, at this scale riccions heavier than M cw can not 
survive, so m ew < M cw . Now eq.(5.1) yields 

PA ™ = 4(ky M ™ [3V * " ln{l + ^ )] (5 ' 2) 

taking m ew = M ew . The second quantization and uncertainty relation 
imply that vacuum has energy density as well as pressure [13]. Ex- 
perimental probes, like la supernova and WMAP [1, 2, 3, 37] suggest 
accelerated expansion of the universe, which requires negative pres- 
sure for the vacuum. So, to have consistency with recent probes, the 
isotropic vacuum pressure is calculated as 
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(2ir)~ 3 /" m ° w 7r 

P^ = -'-i-j o 4rf* = (5.3) 

It yields 

Pa cw /pa ow = cu Acw = —[SV^-lnil+V^)}- 1 = -0.397 ~ -0.4. (5.4a) 

Though in certain models, time- dependence of uo = p/p is also pro- 
posed , but normally it is taken as a constant. So, here also, this ratio 
is considered independent of time. As a result 



= Pa/ Pa = Pa cw /pa ow = -0.4. (5.46) 

Using = M cw = mew and connecting eqs.(4.7),(4.8) and (5.2), a 
quardatic equation for Aew is obtained as 



(450 - Aew 2 + 20 (450 - + 30 x 2818.8] Aew + 11375.3 = 0, 

(5.5) 

which is solved to 

A cw = -0.013352958. (5.6) 
Also using M ew = mew in eq.(5.2), it is obtained that 

p Acw ~ 10 6 GeV 4 . (5.7) 

Planck scale is supposed to be a fundamental scale in field theories. 
So, it is proposed that energy mass scale \x falls from the Planck mass 
M P = 10 19 GeV to the cut-off scale M ew = lOOGeV. When it happens 
so, phase transition takes place at /i — M ew and energy with density 

p cw(r) = A - A cw = 2.5 x 10 7 GeV 4 (5.8a) 



2, 



8tt 2 
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is released, which is obtained connecting eqs.(4.3a) and (5.6). 

In the proposed model (PM), this event is recognized as big — bang, 
being beginning of the universe like standard model(SMU) with the 
release of background radiation having energy density p ew (r), given 
by eq.(5.8a). In contrast to SMU, here, released energy density at 
the epoch of big-bang is finite (it is infinite in SMU). Temperature of 
photons T cw with energy density p cw , is obtained from 

Pew(r) = ^fl = 2.5 x 10 7 GeV 4 (5.86) 

as 

T ew = 78.5GeV = 9.1 x 10 14 K (5.9) 

6. Proposed cosmological scenario, Dark energy and Dark matter 

In the standard model of the big — bang theory, it is supposed that, 
around 13.7Gyrs ago, there used to be a fireball ( an extremely hot 
object), which was termed as primeval atom by Lemaitre. Our uni- 
verse came into existence, when this primeval atom burst out. This 
event is called big — bang. 

According to the proposed cosmological picture, the main content of 
the primeval atom was riccions,being contribution of 10-dimensional 
higher- derivative gravity to the 4-dimensional world. One-loop renor- 
malization of riccions and solutions of resulting group equations yield 
that when energy mass scale comes down to electroweak scale M ew = 
lOOGeV, riccion contributes dark energy density pa ow = 10 6 GeV 4 . 
Moreover, phase transition takes place at M ew , releasing the back- 
ground radiation. This radiation thermalizes the universe upto the 



20 S.K.SRIVASTAVA 

temperature T cw = 9.1 x 10 1A K. Here, the event of phase transition is 
recognized as big — bang, which heralds our dynamical universe having 
the initial temperature T ew = 9.1 x 10 14 -K" and initial value of dark 
energy density p Aow = 10 6 GeV 4 . 

Moreover, it is important to mention that existence of riccions are 
possible at energy scales where higher- derivative terms in the gravita- 
tional action (2.1a) has significant role compared to Einstein-Hilbert 
term. At energy scales below M cw , Einstein-Hilbert term dominates 
higher- derivative terms in the action (2.1a). So, riccions has no di- 
rect role in the evolution of proposed model of the universe, but it has 
two very important contributions to the observable universe as value of 
dark energy density p\ cw = 10 6 GeV 4 and cosmic background radiation 
temperature T ew = 9.1 x 10 14 i^ at cosmic time t — 0. 

Astronomical observations have compelling evidences that the cur- 
rent universe is dominated by dark energy (DE). The present cosmic 
dark energy density is very low, but it used to be very high in the early 
universe. The fall of dark energy density from very high to extremely 
low value can be explained if it is time-dependent. So, like other other 
cosmic dark energy models, here also, dark energy density pa is slowly 
varying function of time. Following Bronstein's idea [14, 15, 16] that 
DE could decay to hot or cold dak matter, here it is proposed that p A (t) 
decays to hot dak matter(HDM) till temperature is high and cold dark 
matter after decoupling of matter from radiation. As a result, p\(t) 
falls from high value 10 6 GeV 4 to currently low value 0.73p cri0 (where 
p CT $ is the critical density). 

It is proposed that the homogeneous dynamical universe begins with 
topology having the distance function 
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dS 2 = dt 2 - a 2 (t)[dx 2 + dy 2 + dz 2 } 



(6.1) 



for the spatially flat model of the universe supported by recent exper- 
iments [1, 2, 3, 34]. This space-time is a special case of hypersurface 
M 4 of the line element (1.1). Here a(t) is the scale factor. 

In what follows, it is obtained that dark energy density decreases 
with time from its initial value 10 6 GeV 4 and falls down by 53 orders 
in the present universe. In 1933, Bronstein proposed that, in the ex- 
panding universe, dark energy density decreases due its decay as a 
result of emission of dark matter or radiation [15, 16]. The radiation, 
so emitted, could disturb spectrum of 3K-microwave background radi- 
ation. So, Bronstein's original idea was modified and it was introduced 
that dark energy (DE) could decay to hot or cold dark matter with- 
out any harm to spectrum of 3K- microwave background radiation [15]. 
Following this idea, here, it is demonstrated that dark energy decays 
to dark matter providing a solution to cosmic coincidence problem. 

It is shown above that , in the beginning, the universe was very hot 
due to the background radiation (released during phase transition). 
Radiation energy density falls as 



with growing scale factor. 

Matter remains in thermal equillibrium with the background radia- 
tion for sufficiently long time. According to WMAP [37], decoupling 
of matter from background radiation takes place at tj ~ 386kyr = 
1.85 x 10 37 GeV _1 . Equation of state for radiation is uo = 1/3. So, here, 



Pr = 



'ew 



(6.2) 



a±(t) 



22 S.K.SRIVASTAVA 

it is proposed that dark energy decays to HDM till decoupling time tj 
obeying uj m = 1/3. 

When t > td, it decays more to CDM ( which is non-baryonic and 
pressureless) with oo m = as well as HDM. So, ratio of densities of 
HDM and CDM is extremely small below td- 

In what follows, development of the universe is probed using these 
ideas. 

(a) Decay of dark energy to dark matter 
The conservation equations 

T (A)j;i + T (dm)j;i = ( 6 -3) 

yield coupled equations 

Pa + 3if (1 + uja) P a = -Q(t) (6.4a) 

and 

p dm + 3H(1 + u dm )p dm = Q(t) (6.46) 

with p\{t) and pd m (t), being dar/c energy density and dark matter 
density respectively at cosmic time t. Here H = a/a,UA = Pa/ Pa = 
—0.4 and uj dm = Pdm/pdm • Q(t) is the loss (gain) term for DE(DM) 
respectively. 

Batchelor [28] has pointed out that dissipation is natural for all ma- 
terial fluids, barring superfluids . So, presence of a dissipative term for 
dark matter DM is reasonable. In eq.(6.4b), this term is obtained by 
taking 
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Q(t) = 3nHp dm (6.5a) 

(with n being the real number) as in [29] without any harm to physics. 
With this setting for time-dependent arbitrary function Q(t), eq.(6.4b)is 
obtained as 

p dm + 3H(1 - n + w)p dm = 0. (6.56) 



Here — np dm given by Q(t) acts as dissipative pressure for DM. As Q(t) 
is proportional to p dm , this setting does not disturb the perfect fluid 
structure as required by the 'cosmological principle'. Thus eq.(6.5b) 
yields the effective pressure for DM as 

P(dm, e ff) = (-n + w)p dm . (6.5c) 

Now connecting eqs.(6.4b) and (6.5) and integrating, it is obtained 
that 

p dm = Aa 3(n - Wdm - 1) (6.6a) 

with 

A = 0.23p cr , o ao 3(n_1) (6.66) 

using current value of dark matter density p dm) o = 0.23p crj0 and a = 
a(t ). 

Connecting eqs.(6.4a), (6.5) and (6.6a) and integrating, it is obtained 
that 

D n 

~ a 3(l+^ A ) _ n + u j A - u± 

where 

D = 10 6 a^+ WA )GeV 4 



(6.7a) 
(6.76) 
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using the initial condition pdm = at t = 0. 
(b) Expansion of the universe, when p A > p dl) 
Friedmann equation is given as 



/Cl\ 2 8hGn / 07TLT7V 
y-J = o (^A + Pdm + Pr) ^ 7, \PA + Pdn 



8nG, 



2 vrJ1 . ruill . r ,, — g vrJ , . ruui/ (6.8) 

as (pA + Pdm) dominates over p r , which is clear from eqs.(6.2),(6.6) and 
(6.7). 



In this case, the equation (6.8) reduces to 



H 2 



_ 8ttG n D 
a) ~ 3a 3 ( 1+ ^A) 



yielding the solution 



a(t) = a c 



1 + 



8nG N D t 
3 

Uew 



10/9 



= a r 



10/9 



1 + 2.89 x l(T 16 t 
with cja = —0.4, a cw = a>(t = 0) and 

= 2.986 x 10- 16 ai( w GeV 



\nG N D 



(6.9) 



(6.10a) 



(6.106) 



( G N = Mp 2 , M P = 10 19 GeV). 

It shows an accelerated expansion of the universe from the beginning 
itself as a/a > 0. 

Connecting eqs.(6.7) and (6.10) 



Pa = 



10 e 



n 



But, the present universe obeys the condition 



"Pdn 



(6.11) 
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aY,o + On,o = l, (6.12a) 

where = Q dmfi + Q rfi . Here Q = p /p C r,o and p crfl = 3H$/8nG N 
(H being Hubble's constant) for the present universe. With these 
values the present critical density is calculated as 

p cr>0 = 3H$/8irG N ~ 1.2 x l(T 47 GeV 4 (6.126) 

using if = ho/ to with /i = 0.68 and the present age of the universe 
to = 13.7Gyr = 6.6 x lO^GeV" 1 [37]. 

As, in the present universe, p m is dominated by cold dark matter 
density, oj m = is taken in eq.(6.11). Now, conecting eqs.(6.11) and 
(6.12), it is obtained that 

n = 0.47 (6.13) 

Using eq.(6.10), the scale factor at the decoupling time t d = 1.85 x 
10 37 GeV _1 and the present scale factor a are obtained as 

a d ~ 1.38 x 10 24 a ew (6.14) 

and 

a ~ 1.58 x 10 29 a ew (6.15) 

using the present age of the universe given above. 

The equation(6.10) shows an accelerated growth of the scale factor 
when pa > pd m exhibiting non-adiabatic expansion of the universe. 
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It implies non-conservation of entropy of the universe, which is in con- 
trast to the decelerated adiabatic expansion of SMU, when cosmological 
model is radiation-dominated or matter-dominated. 

In what follows, scale factor dependence of temperature and en- 
tropy is obtained in an emperical manner, based on current temper- 
ature of the microwave background radiation To = 2.73i^, the ini- 
tial temperature T ew = 78.5GeV = 9.1 x 10 14 -fT given by eq.(5.9) and 
a — 1.58 x 10 34 a cw given by eq.(6.15). Thus , the required emperical 
relation is obtained as 

r -1 103/50 n 

T/T ew ] =*Z. (6.16) 
which yields the decoupling temperature as 

50/103 



T d = T ew — =1740K, (6.17) 

L a d J 

using eqs.(5.9) and (6.14). This value is much lower than = 3000i^, 
obtained in the standard big - bang cosmology. These drastic changes 
are due to dominance of the dark energy in the proposed model. 
Using eq.(5.8b), entropy of the universe is calculated as 

S=^T\ (6.18) 

Current value of entropy S is supposed to be 10 87 . "How could so 
high entropy of the universe be generated ?" is an old question. A so- 
lution to this problem was suggested in a seminal paper on inflationary 
model of the early universe by Guth [38] and modified version of the 
same by Linde [39] and Albrecht and Steinhardt [40]. Here, a differ- 
ent answer to this question is provided on the basis of results obtained 
above. 
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Connecting eq.(6.15) and eq.(6.18) and current vlues of entropy as 
well as temperature T = 2.73-K" of the universe, it is obtained that 



a ew = 0.25 (6.19) 



as 



7vr 2 

S = 10 87 = -^Tj-U- 58 x 2-73 x 10 29 a cw ) 3 

Connecting eqs.(6.16) and (6.18) it is obtained that entropy grows 
with the scale factor a(t) as 



S = ^T c 3 w a^ 103 « 159/103 (6.20a) 

with the initial value 



S cw = 9 x 10 42 (6.206) 
The equation (6.6) yields the rate of production of HDM 

a 1 - 59 

Phdm = -5.9 x 10- 64 ^ (6.21) 

using ujdm = 1/3 for hot dark matter. This equation shows that rate of 
production of HDM increases with growing scale factor. So, production 
of HDM is responsible for increasing entropy in the universe. In SMU, 
number of photons decide entropy of the universe. But, in PM, photons 
and HDM both are responsible for entropy. As energy of HDM increases 
due to its production, owing to decay of DE, entropy is generated in 
this model. It is unlike SMU, where entropy remains conserved. 

Connecting eqs.(6.6) and (6.16), temperature dependence of pdm is 
obtained as 
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3(n-w dm -l) 
U-ew 



0-23p crj0 



T cw /T 



6.18(ra-w dm -l) 



(6.22) 



Putting w dm = 1/3(0) for HDM (CDM), in eq.(6.22), ratio of CDM 
and HDM densities,p c d m and phdm, is obtained as 



2.06 



Pcdm/ Phdm "ew 

T ew /T (6.23) 

Eq.(6.23) shows that creation of HDM is higher at high temperature. 
But as temperature falls down, creation of CDM supersedes the pro- 
duction of HDM ( production of dark matter owing to decay of dark 
energy). For the current universe, the ratio of eq.(6.23) is obtained as 



Pcdm,o/Phdm,o = 2.1 x 10 29 (6.24) 

using numerical values of a ew ,T ew and T (given above) in eq.(6.23). It 
shows that currently, Phdm is almost negligible compared to p c dm- 

The scale factor a S( j , upto which dark energy vanishes due to decay, 
are obtained as 



a S( j = 6.78ao (6.25) 

connecting eqs.(6.6),(6.7) and (6.10) as well as using p A = 0. Here 
Wd m = is taken as for t > t , dark matter content is expected to be 
dominated by CDM. 

The ratio of densities of dark energy, p A , and dark matter, pdm, is 
obtained as 

/ \ - 21 47 
PA/Pdm = I0{a /a) -—. (6.26) 
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It is interesting to note, from eq.(6.26), that the gap between dark 
energy density and dark matter density decreases with the growing 
scale factor. Connecting eqs.(6.6) and (6.7) and using n from eq.(6.13). 
This equation shows that pdm/ Pa < 1 for 



a ew < a < 3.44a and < t < 3.03t - (6.27a, 6) 

The result, given by eq. (6.27b), solves the cosmic coincidence problem 
[24] . This approach, for solution of coincidence problem has an advan- 
tange that no scalar field is required to represent the dark energy, 
rather it is the contribution of higher-dimensional higher-derivative 
gravity to the observable universe. It is unlike the case of work in 
refs. [25-28], where assumed quintessence scalars are used. 

(c) The universe, in case p A < p dm 

Like eq.(6.27), employing the same procedure, it is also possible to 
find that p dm / Pa > 1 for 

a > 3.44a and t > 3Mt . (6.28a, 6) 

So, for t > 3.03to, the the Friedmann equation is written as 

W = ^^(/ ? A+Pdm+Pr)^^^Pdm^0.27/J 2 (^ao/aJ (6.29) 

using eqs.(6.6),(6.8), (6.13) and definition of p cri0 (given above). Eq.(6.29) 
yields the solution 

= [4<dm + 0-43 x 10~ 43 (t - t A <dm)J , (6.30) 

where aA<d m = 3.44a and tA<dm = 3.03to- Eq.(6.30) also shows an 
accelerated growth of the scale factor, even though p A < p dm . It is 
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interesting to see that this expansion is faster than the expansion in 
the interval < t < 3.03to- 

From eq.(6.25), the scale factor upto which dark energy vanishes is 
a s d = 430ao- So, eq.(6.30) yields the corresponding time as 



t sd = 3.7t = 50.69Gyr. (6.31) 

It shows that decay of dark energy will continue in the interval 
3.03t < t < 3.7t also. So, during this time interval, the dark matter 
will follow the rule, given by eq.(6.6). 

Using eqs.(6.16),(6.18), (6.25)and (6.31), T sd and S sd are calculated 

as 

/ an \ 50/103 

T sd = T ( — - ) = 0.395T = 1.08AT (6.32) 
Va sd / 

Dark matter density at t — t sd is obtained as 



Pdm(sd) = 1-32 x l(T 49 GeV 4 (6.33) 

using a and a sd in eq.(6.6). 

As dark energy vanishes when t > t sd , content of the universe will be 
dominated by CDM ,which is pressureless non-baryonic matter obeying 
the conservation equation 

Pdm + 3ifp dm = 0. 
This equation yields scale factor dependence of p dm as 



,3 2.54 xlO 39 

for t > 3.7t . 



Pdm = Pdm(sd) [a sd /a(t) j = ^— (6.34) 
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Beyond the age of the universe 3.7to, the Friedmann equation looks 
like 



a\2 8nG N 212.8 
- Pdm = ^T7r 6.35 

a) 3 a 6 (t) 

using eq.(6.34). 

Eq.(6.35) yields the solution 

T q/9 "1 2 /3 

a(t) = \a% 2 ± U.6(t - t sd )J , (6.36) 

with a s d and t s d given by eqs.(6.25) and (6.31) respectively. 

On taking (+) sign in eq.(6.36), decelerated expansion of the universe 
is obtained beyond t > 3.7to continuing for ever. But the (-) sign, in 
eq.(6.36), exhibits a contracting universe beyond t > 3.7t . So, ulti- 
mately the contracting universe is expected to collapse to a very small 
size with scale factor, possibly equal to a cw = 0.25. From eq.(6.36), 
with (-) sign, the collapse time is calculated as 



3/2 

ted = t sd + ^— = 18.08t = 247.73Gyr (6.37) 
14.6 

using a S( j = 6.78a and t sd = 3.7t . 



7. Elementary particles, Primordial nucleosynthesis and 
Structure formation 

7. 1 Creation of particles 

In SMU, it is assumed that elementary particles such as leptons, 
mesons, nucleons and their anti-particles were produced at the epoch 
of big-bang. Here, production of scalar and spin-1/2 particles is pro- 
posed in the very early universe due to topological changes caused by 
expansion of the universe. 
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A lot of work has been done in the past, where production of particles 
are discussed in curved space-time. This mechanism is based on the 
fact that, in the flat space-time, vacuum state in the Fock space is 
stable, if there is no source term in the action of the matter field. In 
curved space-time, vacuum is unstable even in the absence of the source 
term in the action. It turns out that gravity makes the initial vacuum 
state |0, in > unstable such that |0, in >^ |0,out >, where |0, out > 
is the other vacuum state in the new Fock space. The instability of 
vacuum state shows production of particles where conformal symmetry 
is broken. Moreover, changing gravitational field contributes rest mass 
to produced particles [22, 32, 41]. Here |0, in > state is defined when 
t — > and |0, out > is obtained for t — > oo. 

(a) Creation of spinless particles 

The scalar field obeys the Klein-Gordon equation 



where is mass of 0. Expanding in terms of mode k, eq.(7.1) is 
written as 



(□ + mJ)0 = O, 



(7.1) 



(1 + 2.89 x l(T 16 t)0 fe + 



28.9 



x 1O~ 16 0+ mJ fc (l + 2.89 x 1(T 



3 




'ew 



(7.3) 



using the scale factor a(t), given by eq.(6.10a). 
For small t, eq.(7.3) is obtained as 




(7.4a) 
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where 

r = 1 + 2.89 x l(T 16 t. 
Eq.(7.4) yields the normalized solution 



33 



(7.46) 



where 



and 



= mnf^r^er^F, (| + ^ 2^) 

[ 1+ (zgL+Ml], ( 7 .5a) 
20 (10 15 - 5 A;) 2 



[2(27r) 3 v / 62]" 1/2 e- ir 



61 = 



9 < 



(7.56) 



6 2 = 10 30 (mJ + 



life 2 



Here i-F\(a, 6, c) is the confluent hypergeometric function. 
For large t, eq.(7.3) reduces to 



^%« + (10-WrC'=0, 



dr 2 3 dr 
which integrates to 



(7.5c) 



(7.6) 



out _ r -7/6 



Ci J_ 7/6 (At) + c 2 Y_ 7/6 (At) 



where 



(7.7a) 



A = 10 15 - 5 m^. (7.76) 

and J p (x) and Y p (x) are Bessel's function of first and second kind re- 
spectively. For large x, Bessel's functions can be approximated as 

cos(x — it/A — pn/2) 



J p (x) ~ 



^nx/2 



(7.7c) 
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and 

Y p (x) - gm(x "^lZ W2) . (7.7d) 
Using these approximations, when r is large, eq.(7.7a) looks like 

— 5/3 

0r - -j==[c iC os(At + ^) + c 2 sin(Ar + |)] (7.8) 

Solutions (7.5a) and (7.8) yield number of produced spinless particles 
(for mode k) per unit volume as 

-10/3 

where 

X = [( ~ 3bl ^ T - ^Vh)cos(Ar) - Ay/hTsin(AT)\ +i[{y/b 2 r 

(7.96) 

is defined in Appendix B. Non-zero \f3k\ 2 shows creation of spinless 
particles. 

(b) Creation of spin- 1/2 particles 

The spin-1/2 field ip satisfies the Dirac equation 



(i^D^- m f )^ = 0, (7.10) 

where is mass of ip, 7 M are Dirac matrices in curved space-time and 
are covariant derivatives defind in Appendix B. ip can be written 

as 

*l> = Yl J2( b k,s^ik,s + d\^nk,s) (7.11) 

s=±l k 

with bk yS and c4 )S , given in Appendix B. Eqs.(7.10) and (7.11) yield 
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{i^D^ - m f )ip Ik)9 = 0, 



(7.12a) 



(iYD li -m f )^ IIk , 8 = 0. (7.126) 

Now using the operator (—ij^D^ — rrif) from left of eqs.(7.12a) and 
(7.12b), it is obtained that 



;□ + -R + m 2 f )j> = 0, 



where ip = ^i k , s {^iik,s)- 
Writing 

and 

eq.(7.13) is obtained as 



ipik, s = fik, s (tV k '* u 



(7.13) 



d 2 f 10 df 



+ 10 



30 



k 2 4.93 x 10" 15 



m f T — 



f T ll/9 



dT 2 3 dT 

where /= fi k , s (fnk,s)- 

For small t , eq.(7.14) looks like 



3r 



f = 0, (7.14) 



d 2 f 10 df in30 



f 2 _|_ ■ 



T — 



20k 2 

Q/,2 



7=0. (7.15) 



This equation is like the equation (7.4a). So, its solution has the 
same form. Using this solution in tpzk,s defined above , it is obtained 
that 



As = 



1 H — ' 



10 



(7.16) 



For large t , eq.(7.14) reduces to 
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yielding 



_d 2 f 10 df 

d^ + Yd^ 



+ 10 30 mjr/ = 



(7.17) 



-5/3 



IIk,s 



7T 2 A 



■cos(Ar)e lk - x u s . 



(7.18) 



Using i^Yk,s an d i'Tik^ m ^ ne definition of /3 fe;S , given in the Appendix 



B. 



„ ^3^-5/3 ^^) t 



(_ 3&1 + 10 ?v ^) 

1 + To r 



(7.19) 

So, the number of created spin-1/2 particles per unit volume is ob- 
tained from eq.(7.19)as 



\Pk,s\ 2 = 



(27T) 2 



-10/3 



3&1 V 1 h 2 

T + OoT 

10 ' 



(7.20) 



2 7 T 2 ,4v / &2 

Eqs.(7.9a,b) and (7.20) show creation of spinless and spin-1/2 par- 
ticles due to changing gravitational field in the proposed speeded-up 
model. The scale factor, given by eq.(6.10a) show that the significant 
change in the gravitational field is possible , in this model, when 



* > 3.46 x 10 15 - 5 GeV _1 = 7.2 x 10" 9 sec. (7.21) 

as a(t) remains almost constant upto this epoch. So, here, particle 
production is expected, when the universe is around 2.3 x 10" 9 sec. old. 
Also, these results show that particle production falls down rapidly as 
time increases. Thus, like other models, here also particle production 
is expected in the early stages of the universe. 

7.2 Primordial Nucleosynthesis 
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Hydrogen (H) is the major component of baryonic matter in the 
universe. The next main component is Helium-4 { 4 He) . Occurrence 
of other light elements and metals is very small. It is found unlikely 
that abundance of A He, Deuterium (D) and other light elements, in 
the universe, being caused by burn out of H in stars [42, 43] So, it is 
argued that the required amount of these might have produced in the 
early universe. 

At the epoch of its formation, helium production depends upon neu- 
tron (n) concentration, which is determined by weak interaction reac- 
tions given as 



n+u^p+e ,n + e + ^p + u (7.22) 

where p stands for proton and v for neutrino. This chemical equilib- 
rium is maintained till weak reaction rate T w >> H, where H is the 
expansion rate of the universe given by eq.(6.9) and T w ~ 1.3G 2 F T 5 
with Fermi constant G F = na w /V2Ml = 1.17 x l(T 5 GeV~ 2 . With the 
expansion of the universe, temperature decreases so weak interaction 
rate slows down . As a result, at the freeze-out temperature T*, 

T w ~ H. (7.23a) 

Connecting eqs.(6.7b), (6.9), (6.16) and (7.23a), it is obtained that 

1.3 x (1.17) 2 x KT 10 ?; 5 = 2.89 x 1(T 16 (^-) ^ 
which yields 

T* = 0.9MeV. (7.236) 

Eqs.(6.10a,b) and (6.16) yield time dependence of temperature , for 
this model, as 
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4.5 

^sec = rpi 9 • (7.24) 
J MeV 

Using this result, the freeze-out time t*, corresponding to T*, is ob- 
tained as 

U ^ 5.55sec. (7.25) 

In SMU,T, ~ 0.86MeV and t* - ^sec [44]. It is obtained that, in 
PM, T* is a bit higher and freeze-out takes place later. The reason 
for these differences is the basic difference in developement of these 
models. SMU is driven by elementary particles present in the early 
universe and it expands adiabatically. The proposed model is driven 
by dark energy, with very high density at the beginning and it expands 
with acceleration. 

Like [45], in this subsection, temperature T 9 is measured in units 
10 9 fT, so lMeV becomes 11.6 in T 9 temperature. Now, eq.(7.24) looks 
as 



474 

* = ^U- ( 7 - 26 ) 
For large temperature T 9 > 10, neutron abundance is given by 



X„ 



1 +e 



Q9/T9 



-1 



(7.27a) 



where Q 9 = m n — m p = 15. At freeze-out temperature T* ~ 0.9MeV, 
which is equivalent to T 9 = 10.44, 



X* = 0.19. (7.276) 

Due to slightly higher freeze-out temperature, in the present model, 
neutron concentration is obtained higher than X* = 0.16 in SMU. 
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When r w << H i.e. when chemical equilibrium between n and p 
frezze-out, neutron concentration is determined through neutron-decay 

n^p + e~ + v. (7.28) 

Thus, for t > £*, neutron abundance is given as 

X n (t) = 0.19e^ /r , (7.29) 

where r = (885.7 ± 0.8)sec. is the neutron life-time. 
First step, in the formation of complex nuclei, is 

p + n <-> D + 7 

showing chemical equilibrium between deutron (D) , nucleons and 
photon (7). This equilibrium gives deutron abundance by weight [43] 
as 

X D = 0.2 x lO- 12 (n B h 2 )X n X p T* /2 exp(B D /T 9 ), (7.30) 

where X D = 2n D /n B ,X n = n n /n B ,X p = n p /n B ,X n + X p — 1, deu- 
terium binding energy Bp = m p + m n — mo = 2.23MeV which yields 
Bog = 25.82 (in units of Tg) and (Q B h 2 ) is the baryon number density. 
Light heavy elements 4 He is formed through reactions 

n + p^D + 'y, (7.31a) 

(i)D + D ^ 3 He + n, (ii)D + D — > 3 H + p (7.316, c) 

and 

(ii)D + 3 H — » 4 He + n, {Of He + D^ 4 He + p, (7.31d, e) 
where 3 if stands for tritium. 
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These reactions show that sufficient deuterium abundance is required 
for nucleosynthesis of 4 He, which is supposed to be around 25% of the 
baryonic content of the universe. But after weak-interaction freeze out, 
deuterium abundance is not sufficient unless temperature is very low. 
For example, at T 9 = 5.8(T = 0.5MeV),X D = 2 x KT 12 is obtained 
using eq.(7.30) and putting fish 2 = 0.05. 

Rates for reactions (7.30b, c) are given as [45, Appendix] 

\ DD(i) = 3.97 x 10 5 T 9 " 2/ V 4 - 258/T 9 1/3 [1 + 0.098T 9 1/3 + 0.876T 9 2 /3 

+0.6T 9 - 0.041T 9 4/3 - 0.071T 9 5/3 ](fi i? /i 2 )sec- 1 
Xdd(u) = 4.17 x 10 5 T 9 _2/ V 4 - 258/T 9 /3 [1 + 0.098T 9 1/3 + 0.518T 2/3 
-0.355T 9 - 0.010T 9 4/3 - 0.018T 9 /3 ](fi B /i 2 )se C - 1 

(7.32a, b) 

Making approximations in eqs.(7.30a,b), it is obtained that 

Add = ^DD(i) + ^dd(m) 

~ 8.14 x 10 5 T 9 - 2/3 e- 4 - 258 / T 9 1/3 (fiD/i 2 )sec- 1 

(7.33) 

Mukhanov [44] has obtained the condition for conversion of sufficient 
amount of D to 3 He and 3 H as 

1 

2 

From eqs.(7.27), (7.30)-(7.33), it is obtained that 



2 \ DD X D t~l (7.34) 



25.82 4.258 



6.43 x 10" 6 T 9 7 L3 e V T9i (Toi) 1 ' 3 J (Q B h 2 ) 2 ~ 1 (7.35) 

showing a relation between Q B h 2 and T 9i . T 9i is the temperature at 
which nucleosynthesis takes place. For VL B h 2 = 0.05, this equation 
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yields 

T 9i ~ 1.16, (7.36) 

which corresponds to Tj = O.lMeV. It shows that if baryon number 
density Vtsh 2 is ~ 0.05, nucleosynthesis of 4 He may take place at tem- 
perature Tj = O.lMeV. In SMU, this temperature is ~ 0.086MeV. 
Eqs.(7.29) and (7.35) yield deuterium concencentration at T; as 

X Di ~ 8.93 x 10~ 6 (7.37) 

It happens at the cosmic time 

ti = 357.5sec. (7.38) 

Thus, it is obtained that, in the present model, nucleosynthesis of 4 He 
begins much later compared to SMU, where it happens at time ~ 
lOOsec [44]. But, it is still less than life-time of a neutron, which is 
required for the necleosynthesis. 

The final 4 He abundance is determined by available free neutrons at 
cosmic time tj. As total weight of A He is due to neutron and proton, 
so its final abundance by weight is given by [45] as 

X{ = 2X* n exp(-tf /t) = 2X>p(-0.535/(T/) 1 - 9 (7.39) 

using eqs.(7.26) and r = 886.5sec. Here a stands for 4 He. 
The rate of production of 4 ife is given by the equation 

X a = 2\ DD ^Xfy. (7.40) 

As free neutrons are captured in 4 He, helium production dominates 
neutron decay. It happens when [45] 



2X a = ^. 
r 



(7.41) 
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Approximating Xdd(u), given by eq.(7.32b), around T 9 = 1.16, it is 
obtained that 



\ DD(ll) ~ 6.98 x 10 5 T 9 - 2/3 e^(-4.258/T 9 1/3 )(fi B / i 2 ). (7.42) 
Eqs.(7.29), (7.38)-(7.40) yield 



| (S^Si-l^) , 0.535 v 

— ~ 6.98 x 10 T q 1 e T s x 0.19e T 9 
2r 9 

x(0.2 x 0.81 x 10- 12 ) 2 (n B h 2 ) 3 , 

which implies that 



T f - ^ (7 43) 

J 9 37.96 _ UnT 9 - 3ln(Q B h*) ■ I #§r ' 1 J 

At T = T 9i = 1.16, 



T f = (744 ) 

9 42.08 - 3ln(Q B h 2 Y 1 ' ; 

Using eq.(7.44) in eq.(7.39), helium abundance by weight, as function 
of baryon number density (Q B h 2 ),is obtained as 

t r f 42.08 -3ln(n B h 2 )-) 1 -v 

X{ ~0.38exp[- 0.535 1 ^ J j .. (7.45) 

For Vl B h 2 = 0.05, it is obtained that 

X{ ~ 0.23. (7.46) 

For T < T/, neutron abundance is governed by X n = -2X a with 
X a given by eq.(7.39). So, 

A n = — 2A/)D(m)A^, 
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where Xo and \dd(u) are given by eqs.(7.30) and (7.41) respectively. 
Using eq.(7.26) and integrating, this equation yields neutron abundance 
below T( as 



X„ 



(1/X") + (Q B h 2 ) 3 e- 37 - e2 T° A3 (e 51 - 64/T9 - e 51 - 64 ^)] \ (7.47a) 



which shows that concentration of free neutrons drops very fast. It is 
due to capture of free neutrons in 4 He. From eq.(7.29), it is obtained 
that neutron and deuterium concentrations become equal at tempera- 
ture Tg given as [45] 

T d = 25,82 (7 A7h) 

9 29.45 -ln{Sl B W)-\lnTf K ] 

Deutron concentration is given by [45] 

X D ~ -2\ DD(U) X 2 D = -2 x 8.14 x ltf T - 2 l\n B h 2 )e-^' T » * X* , 
which integrates to 

X D = [(l/X D (T 9 rf )) + (^a^e 20 - 16 !^)- 2 - 5 ^" 4 - 258 /^) 173 
-(T 9 )- 2 - 8 e- 4 - 258 /( T9 ) 1/3 }] _1 , 

(7.47c) 

where X D (T 9 d ) = X n (T$). 

Eqs. (7.47a) and (7.47c) show that, below T/, abundance of free neu- 
tron and deuterium drop rapidly in PM also like SMU. 

7.3 Growth of inhomogeneities in the proposed model 

For small inhomogeneities, linear perturbation of Einstein equations 
is enough. In this case, contrast density 8 = opdm/Pdm (with pdm, being 
the cold dark matter energy density and opdm is small fluctuation in 



44 S.K.SRIVASTAVA 

Pdm presenting inhomogeneity) , obeys a linear differential equation in 
the homogeneous model of the universe [46]. Due to linearity, S can be 
expanded in modes k as 



(7.48) 



For modes k with proper wavelength F < H 1 (t)(H being the 
Hubble's radius), the perturbation equation looks like 

'k 2 v 2 



d fc + i-Ok + 



(-^-47rGp dm )<J fc = 0, (7.49) 



a v a z 

which is free from gauge ambiguities [46, eq.(4.158)]. So, Sk depends 
on t only. 

Eqs.(6.5c) and (6.13) give effective pressure for DM as 



It yields 



P(dm,eff) = (-0.47 + w dm )p dm . 



2 ^P(dm,eff) „ 

f „ = — ; = —0.47 + w. 



s 



dp, 



dm 



(7.50) 



(7.51) 



Connecting eqs.(6.6), (6.13), (7.49) and (7.51), the differential equa- 
tion for Sk is obtained as 



k 2 



(-0.47 + w dm )^- - 0.345# 2 4 5 7a 3 (°- 53+WdE 



5 k + 2-5 k + 



with p crj0 = 3H 2 /8nG. 

Using eq.(6.10), eq.(7.52) looks like 



s k = o 

(7.52) 



^(P6k Wd6k 1q31 



dr 2 3 dr 



0.345# 2 < 59 



-0.47 + w dm )- 



A' 2 



r 10(0.23+w dm )/3_ 

where r is defined in eq.(7.4b). 



r ll/9 

5k = 0, egno(7.53) 
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As mentioned above, for CDM Wd m = 0, so eq.(7.52) is approximated 



as 

d 2 5 k lOdh _ 3.45 x 10 29 ff 2 4 59 
T ~d^ + T ~d^ r 5 ^ 
which is integrated to 



6k = 0, (7.54) 



10/21 V*' 1 " ) ^ ^l ±x -110/21 

ci + c 2 ) J 110 /2i(2*r 7 / 60 ) + i(ci - c 2 )r 110 /2i(2^r 7 / 60 ) 



(7.55) 

where i = v^T and 6 2 = 3.45 x 10 30 # 2 aJ- 59 . Here H^\x) and iz£ 2) (a;) 
are Hankel's functions as well as J p (x) and Y p (x) are Bessel's functions. 
Moreover c 1 and c 2 are integration constants. 

As discussed in section 6, CDM dominates HDM in the late universe, 
so structure formation is expected for large r. Using approximations 
(7.7c,d)for Bessel functions, for large r in eq.(7.55)and doing some 
manipulations, it is obtained that 

4^r- a67 e^ /6 °, (7.56) 
which shows growth of structure formation as cosmic time t increases. 

8. Summary of Results 

In contrast to the original big — bang theory, the proposed cosmol- 
ogy answers many basic questions as (i)"What is the fireball?", (ii) 
"How does it burst out?", (iii) "What is the background radiation?" 
and "What are initial values of temperature and energy density?" 
Moreover, the proposed cosmological model is free from initial sin- 
gularity. It is found that the dark energy violates the strong energy 
condition showing "bounce" of the universe, which is consistent with 
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singularity — free model of cosmology [20, 24, 25]. The initial scale 
factor is computed to be a ew = 0.25. It is in contrast to the standard 
model of cosmology SMU, which encounters with singularity having 
zero scale factor, infinite energy density and infinite temperature. 

The present value of dark energy density is supposed to be ~ 7.3 x 
10 _48 Gev 4 , which is 53 orders below its initial values 10 6 Gev 4 . The 
question "How does dark energy falls by 53 orders in the current uni- 
verse?" is answered by the result, in section 6(a), showing that dark 
energy decays to dark matter, obeying the rule, given by eq.(6.7). 
Upto the decoupling time td — 386kyr, matter remains in thermal 
equilibrium with the background radiation, so produced dark matter 
upto td is supposed to be HDM. But when t > td, production of CDM 
is more than HDM. The current ratio of HDM density and CDM den- 
sity is found to be 5 x 10~ 30 . In section 6, it is found that creation of 
HDM raises entropy of the universe upto ~ 10 87 . During the phase of 
accelerated expansion, temperature falls as a(t) -50 / 103 . 

It is found that one of the two types of dynamical changes of the 
universe are possible, beyond 3.7to = 50.69Gyrs, (i) decelerated expan- 
sion and (ii) contraction. If future universe expands with deceleration 
, it will expand for ever. But in the case of contraction, it will collapse 
by 247.73Gyrs. 

As dark energy dominates over dark matter, from the beginning 
upto the time 3.03to, no cosmic coincidence problem arises in the pro- 
posed scenario. Moreover, in the preceding section, some other basic 
problems such as creation of particles in the early universe, primordial 
nucleosynthesis and structure formation in the late universe. In the 
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proposed model, it is shown that particles are created due to topologi- 
cal changes which is unlike the SMU, where it is assumed that elemen- 
tary particles are created at the time of big-bang. It is also shown, in 
section 7.2, that process of primordial nucleosynthesis goes very well in 
the proposed model predicting ~ 23% Helium-4 abundance by weight. 
In section 7.3, it is shown that ,in the late universe, inhomogeneities in 
CDM grow exponentially with time. 

Thus the proposed model is able to provide possible solutions to 
many cosmological problems with prediction for the future universe. 

Appendix A 
Riccion and Graviton 

From the action 



J d 4 xd D yy/-g {4 +D) 



M^R (4+D) 



16tt 



+ tt(4+D)-R^4+Z)) + 



3 6QD + 3) 2 

7(4+D)l-K( 4 +D) -2) ( D+4 ) K (D+4)) 

the gravitational equations are obtained as 



(Al) 



r(2) 



M( 2+D ) 1 (1) 

j_g {Rmn - ^9mn R{4+d)) + ®.(a + d)H MN +~i( i+D )H MN = 0, (A2a) 

where 



H-MN — 2-R;MAT — ^9MN^(4+D)R(4+D) ~ -9MnR(a+D) + 2R(4+D)RmN, 

(A.2b) 

and 

H MN = 3R^MN -^9MN D (4+D)R 2 ^+D) - ^ D _ {-^9MN n (4+D)RlA + D) 
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+2E3(4+d)R(d+4:)Rmn + R 2 ,mn) ~ 2^ MN ^i 4 + D ) ^R^+d^mn- {A. 2c) 

Taking g MN = r] MN + h MN with r/ MN being (4 + D) -dimensional 
Minkowskian metric tensor components and Hmn as small fluctuations, 
the equation for graviton are obtained as 

a (4+D)h M N = (A. 3) 

neglecting higher-orders of h. 

On compactification of M 4 <g> S D to M 4 , eq.(A.3) reduces to the 
equation for 4-dimensional graviton as 

□V + = (A4) 

for the space time 

dS 2 = g^dx^dx" - p 2 [d6l + sin 2 e x de\ H h sm 2 6>i • • • sin 2 9^ D -i)d9 2 D . 

(A5) 

The 4-dimensional graviton equation (A. 4) is like usual 4-dimensional 
graviton equation ( the equation derived from 4-dimensional action) 
only for I — 0. Thus the massless graviton is obtained for / = only. 

As explained, in section 2, the trace of equations (A. 2) leads to the 
riccion equation 

[□ + ±£R + m\ + ^R 2 ]R + -d = 0, (A6a) 

where 
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(D + 2)\V D DRp 1 2 2 
m « " " 16ttG (4+d) + 2(DTZj + 2 V XRd 



A = 



1 



A(D + 3)a 



* ~ 16^ + 4(5T3) + e 77 XRd 

The graviton /i^ has 5 degrees of freedom (2 spin-2 graviton, 2 spin- 
1 gravi- vector (gravi-photon) and 1 scalar). The scalar mode / satisfies 
the equation 

P 

Comparison of eqs.(A.6) and (A. 7) show many differences between 
scalar mode / of graviton and the riccion (R) e.g. £, A and d , given by 
eqs.(A.6b,c,d,e), are vanishing for /, but non- vanishing for R. Eq.(A.7) 
shows {mass) 2 for / as 

2 l(l + D-\). (Ao . 
m) = if, (A8) 

whereas (mass) 2 for R, given by eq.(A.6c), depends on G^ + d,Vd and 
Rd (given in section 2). mj = for / = 0, but m|, can vanish only 
when gravity is probed upto ~ 1CT 33 cm. As mentioned above, so far, 
gravity is probed only upto 1cm. Thus (mass) 2 of riccion does not 
vanish. 

So, even though, / and R are scalars arising from gravity, both are 
different. Riccion can not arise without higher-derivative curvature 
terms in the gravitational action, but graviton can be obtained even 
from Einstein-Hilbert action. 
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1 9 ( , „„ d \ .... d 



^cVA v yy dx v ) ' dX^dX^ 
where X^ are locally inertial co-ordinates and rf v are Minkowskian 

metric components. It shows that the scalar like operator □ has the 

same role on R as it is for other scalar fields due to principle of 

equivalence. 

Appendix B 
Bogoliubov transformations 

In a Hilbert space, a sacalar <3> satisfying the Klein-Gordon equation 
can be written as linear combinations 



n 

k 

J2 [^r$r + 4 out $r ut ] > 



(5.1a) 

where $^ ut and $ k n both belong to the same Hilbert space. As a result, 
one obtains 

where a k and fi km are Bogolubov coefficients satisfying the condition 

K| 2 -|/3 fc | 2 = i. (5.2) 
The in- and out- vacuum states are defined as 

4 n |in >= = A° ut |out > . (5.3a, 6) 

Moreover, 

AT = a k A™ + (3 k Af\ (5.3c) 
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The normalization condition for $ is given as 

($ fe , $ fe ) = 1 = -($*, $*), ($ fe , $*) = (5.4a, 6, c) 

where the scalar product is defined as 



where E is the 3-dim. hypersurface. 
Connecting eqs.(B.lb) and (B.4a,b,c), 

and 

A = -($r,<fD (5.56) 

(b) Spin-1/2 field 

The spin-1/2 field ip, satisfies the Dirac equation 

(i-fD li -m f )il) = ti, (B.Qa) 

where m/ is mass of ^ and 

D, = <9 M - T M (B.6b) 

with Dirac matrices in curved space-time 

7 M = e^7 a , (B.7a) 

where (a, /x = 0, 1, 2, 3) and are defined through 

= Vab- (B.7b) 

Here r/ a fc are Minkowskian metric tensor components and g^ v are met- 
ric tensor components in curved space-time. c.c. stands for complex 
conjugation. 
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Dirac matrices 7 M in curved space-time satisfy the anti-commutation 
rule [22, 32] 

{y, v} = 2*r (b.7c) 

and Dirac matrices 7 a in Minkowskian space-time satisfy the anti- 
commutation rule 

{7 a ,7 fe } = 2 V ab . (B.7d) 

are defined as 

r M = -\ (drf + r;e:)^ff. (s.? e ) 

Further, *0 can be decomposed as [22] 

s=±l fc 
s=±l fc 

(5.8a, 6) 

as both in- and out-spinors belong to the same Hilbert space. The in- 
and out-vacuum states are defined as 

fc£ s |m >= cT M |m >= (5.9a, 6) 

and 

6g|out >= cP^Jout >= (5.9c, d) 
Bogoliubov transformations are given as [22] 



TOUt 

°k,s 


= V£ ta a k>a + d j X- s Pk,s 


i tout 
U k,s 


= <*ljZ + KAk,-. 


/tout 

'-fe,-s 


= K,s a k,s + d^l_ s p k , s 


rout 

'-fc,-s 


Jj^ n _|_ Q* ^ n 



(5.100,6, c, d) 



and 
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Connecting eqs.(B.8)-(B.10), one obtains [33] 

Kl 2 + l^l 2 = l afc -l 2 + I&--I 2 = x > ( 5 - n ) 

s 

where 

a k . s = J^^^x^l s) ^k,s) (5.12a) 

a.- = | s « 3 *; M )7°« M ( 5 - 12fe ) 

Acknowledgements 

Author is thankful to Inter University Centre for Astronomy and 
Astrophysics, Pune for hospitality during my visit in January 2004, 
where idea for this work came to my mind during 3 days seminar on 
brane-world. Thanks are also due to Varun Sahani for helpful sug- 
gestions. Useful discussions with Prof. K.P.Sinha, about riccions, is 
acknowledged gratefully. 



References 

[1] A. Ricss et al., Astron. J. 116, 1009 (1998); astro-ph/9805201; S. J. Pcrlmuttcr 
et al, Astrophys. J. 517, 565 (1999); astro-ph/9812133; J.L. Tonry et al., 
astro-ph/0305008. 

[2] D. N. Spergel et al, Astrophys J. Suppl. 148 175(2003)[ astro-ph/0302209]. 
[3] J. A. Peacock et al., Nature, 410, 169 (2001). 

[4] F. Hoyle and J.V. Narlikar, Proc. R. Soc. London A 282, 191 (1964); Mon. 
Not. R. Astron. Soc. 155, 305 (1972); 155, 323 (1972). 

[5] B. Ratra and P. J. E. Peebles, Phys.Rev. D 37, 3406(1988); C. Wetterich, 
Nucl. Phys. B 302, 668 (1988); J. Frieman, C.T. Hill, A. Stebbins and I. Waga, 
Phys. Rev. Lett. 75, 2077 (1995); P. G. Fcrrcira and M. Joyce, Phys.Rev. D 
58, 023503(1998); I. Zlatev, L. Wang and P. J. Stcinhardt, Phys. Rev. Lett. 
82, 896 (1999); P. Brax and J. Martin, Phys.Rev. D 61, 103502(2000); L. A. 
Urena-Lopez and T. Matos Phys.Rev. D 62, 081302(R) (2000); T. Barrriro, E. 
J. Copeland and N.J. Nunes, Phys.Rev. D 61, 127301 (2000); A. Albrecht and 
C. Skordis, Phys. Rev. Lett. D 84, 2076 (2000); V. B. Johri, Phys.Rev. D 63, 
103504 (2001); J. P. Kneller and L. E. Strigari, astro-ph/0302167; F. Rossati, 
hep-ph/0302159; V. Sahni, M. Sami and T. Souradcep, Phys. Rev. D 65, 
023518 (2002); M. Sami, N. Dadhich and Tetsuya Shiromizu, hep-th/0304187. 



54 



S.K.SRIVASTAVA 



[6] C. Armendariz-Picon, T. Damour and V. Mukhanov, Phys. Lett. B 458, 209 
(1999); T. Chiba, T. Okabc and M. Yamaguchi, Phys. Rev. D 62, 023511 
(2000) 

[7] A. Sen, J. High Energy Phys. 04, 048 (2002); 07, 065 (2002); Mod. Phys. Lett. 
A 17, 1799 (2002). 

[8] M. R. Garousi, Nucl. Phys. B 584, 284 (2000);J. High Energy Phys. 04, 027 
(2003). 

[9] E.A. Bergshoeff et al, J. High Energy Phys. 05, 009 (2000). 
[10] G. W. Gibbons, Phys. Lett. B 537, 1 (2002). 

[11] Parampreet Singh, M. Sami and Naresh Dadhich,Phys.Rev. D 68, 023522 

(2003) and references therein. 
[12] A. Mazumdar, S. Panda and A. Perez-Lorenzana, Nucl. Phys. B 584, 284 

(2001). 

[13] V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. D 9, 373 (2000); astro- 
ph/9904398. 

[14] V. Sahni, Class. Quantum Grav. 19, 3435 (2002); astro-ph/0202076 ; 

S.K.Srivastava, astro-ph/0407048. 
[15] B. Ratra and P. J. E. Peebles, Rev. Mod.Phys. 75, 559(2003), astro- 

ph/0207347. 

[16] M. Bronstcin, Phys. Z. Sowjetunion 3, 73 (1933). 
[17] R. Maartens, gr-qc/0312059. 

[18] A.A.Starobinsky ,Phys.Lctt., 91 B, 99 (1980); Phys.Lett., 157 B,361 (1985) . 

[19] Ya.B.Zeldovich & I.D.Novikov, Rclyativisteviakaya Astrofizika (Relativis- 
tic Astrophysics), Fizmatgiz (1968) ; Comm.Math. Phys., 27, 185 (1971); 
V.Ts.Gurovich & A.A.Starobinsky, Zh.Eksp.Teor.Fiz., 73, 369 (1977) 
[Sov.Physa.JETP,46,193 (1977)] ; A.A.Starobinsky ,Phys.Lett.B, 91 99 (1980) 
; A.A.Kofmann, A.D.Linde & A.A.Starobinsky, Phys.Lett.B, 157,361 (1985) 
; B.Whitt, Phys.Lett.B, 145,176 (1984); S.W.Hawking & J.C.Luttrell , 
Nucl.Phys.B, 247,250 (1984) ; N.H.Barth & S.M.Christensen, Phys.Rev.D, 
28,1876 (1983) ; M.B.Mijic, M.S.Morris & Wai-Mo Sucn, Phys.Rev.D, 34,2934 
(1986) ; Wai-Mo Sucn & P.R. Anderson, Phys.Rev.D, 35,2940 (1987). 

[20] S.K.Srivastava and K.P.Sinha, Phys. Lett. B 307,40(1993) . 

[21] K.P.Sinha and S.K.Srivastava , Jour. Ind. Math. Soc.61,80 (1994). 

[22] S.K.Srivastava and K.P.Sinha , Int. J. Theo. Phys. 35, 135 (1996). 

[23] S.K.Srivastava , Mod. Phys. Lett. A 12, 2933 (1997). 

[24] S.K.Srivastava ,IL Nuovo Cimento,113, 1239 (1998). 

[25] S.K.Srivastava , Int. J. Mod. Phys. A 14, 875 (1999); Int. J. Mod. Phys. A 
15, 2917 (2000). 

[26] P. J. Steihardt, Cosmological challenges for the 21st centuryin : V.L. Fitch, 
D.R. Marlow (Eds.), Critical Problems in Physics, Princeton Univ. Press, 
Princeton, NJ, 1997. 

[27] L. Amendola, Phys.Rev.D, 62, 043511 (2000) ; L. Amendola, D. Tocchini- 

Valcntini, Phys.Rev.D, 64, 043509 (2001) 
[28] G.K.Batchelor, An Introduction to Fluid Dynamics (Cambridge university 

Press, Cambridge, 1967). 
[29] W. Zimdahl, D. Pavon and L. P. Chimcnto, Phys.Lett.B, 521, 133 (2001); 

R. Herrera, D. Pavdn and W. Zimdahl, astro-ph/0404086, Gcn.Rel.Grav. 36 



DARK ENERGY AND DARK MATTER. 



55 



, 2161 (2004); B. Gumjudpai, T.Naskar, M.Sami and S. Tsujikawa, hcp- 
th/0502191. 

[30] L. P. Chimcnto, A. S. Jakubi, D. Pavon and W. Zimdahl, Phys.Rcv.D, 67, 
083513 (2003). 

[31] L. P. Chimcnto and A. S. Jakubi, Phys.Rcv.D, 67, 087302 (2003); Phys.Rcv.D, 
69, 083511 (2004). 

[32] N.D.Birrel and P.C.W.Davies, Quantum Fields in Curved Space (Cambridge 

Univ. Press, Cambridge, 1982). 
[33] S. K.Sri vast ava and K.P.Sinha , Aspects of Gravitational Interactions (Nova 

Science Publications, NY, 1998). 
[34] R.B.Mann, L.Tarasov, D.G.C.Mckeon and T.Stcllc , Nucl.Phys. B 311 , 

630(1988/89). 

[35] B.L.Hu and D.J. O'Connor, Phys.Rev. D 30, 743(1984) . 

[36] E.Elizalde, K.Kirsten and S.D.Odintsov, Phys.Rev. D 50, 5137(1994) . 

[37] A.B. Lahnas, N.E. Mavromatos and D.V. Nanopoulos, Int. J. Mod. Phys. D, 

12(9), 1529 (2003). 
[38] A. H. Guth, Phys.Rev. D 23, 347(1981) . 
[39] A. D. Linde, Phys. Lett. B 114 , 431(1982) . 

[40] A. Albrecht and P. J. Stcinhardt, Phys. Rev. Lett. 48 , 1220(1982) . 
[41] E. Mottola, Phys. Rev.D 31 , 754(1985). 

[42] S. Weinberg, Gravitation and Cosmology (Wiley, New York, 1972). 
[43] E.Kolb and M. Turner , The Early Universe (Addison - Wesley, 1993). 
[44] V. Mukhanov, Int. J. Theo. Phys. 43,669 (2004): astro-ph/0303073. 
[45] R. Esmailzadeh, G.D. Starkman and S.Dimopoulous, Astrophys. J. 378, 504 
(1991). 

[46] T. Padmanabhan, Structure formation in the Universe (Cambridge University 
Press, 1993). 



